Localized edge vibrations and edge reconstruction by Joule heating in graphene 

nanostructures 
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Control of the edge topology of graphene nanostructures is critical to graphene-based electronics. 
A means of producing atomically smooth zigzag edges using electronic current has recently been 
demonstrated in experiments [Jia et al., Science 323, 1701 (2009)]. We develop a microscopic 
theory for current-induced edge reconstruction using density functional theory. Our calculations 
provide evidence for localized vibrations at edge-interfaces involving unpassivated armchair edges. 
We demonstrate that these vibrations couple to the current, estimate their excitation by Joule 
heating, and argue that they are the likely cause of the reconstructions observed in the experiments. 
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Graphene, a single sheet of carbon atoms in a hexago- 
nal lattice, is a material currently under intense scrutiny 
[H, 0|- Graphene is interesting not only because of its 
exotic material properties, but even more so due to its 
potential use in future electronic components. Graphene 
electronics has a tremendous potential 0, but its 
practical realization requires the ability to manufacture 
graphene nanostructures in a controlled and efficient 
manner. The topology of graphene edges plays a funda- 
mental role in determiningthe electronic and transport 
properties of these devices [1,0,0] • Thus, the control and 
stability of edges is crucial for further development of 
graphene-based electronic devices. Recent experiments 
[3, Lol] snow that simple armchair and especially zigzag 
edges are the most commonly occurring edge structures, 
and that their formation and dynamics is strongly influ- 
enced by the energetic electrons in a transmission elec- 
tron microscope (TEM). However, also intermediate re- 
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experiment Jia et al. [12| demonstrated the formation of 
smooth zigzag edges from disordered edges in graphene 
in the presence of an electronic current. The possibility 
of an in situ fabrication process, as suggested by this ex- 
periment, is very attractive. However, at the same time 
the devices should remain stable in the presence of elec- 
trical current for reliable operation, further underlining 
the importance of understanding the microscopic edge 
reconstruction mechanisms. 

In this Letter we present an ab inito study of current- 
induced edge reconstructions in systems, where armchair 
and zigzag edges are adjacent. As the first example, con- 
sider the graphene nanoribbon (GNR) junction shown 
in Fig. [TJ A zigzag-GNR to the left is connected to 
a wider zigzag-GNR to the right by an armchair edge 
(ZAZ-structure). Note that the edges are not passivated 
(we return to this point below) . The second example is a 
ZAZZZ-system (see Fig. [2]) where an extra zigzag edge is 
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FIG. 1: (Color online) [Top] The ZAZ structure. [Middle] 
The two modes lying outside the bulk bands, shown in red 
(a) and blue (b), respectively. Atomic displacements of less 
than 5% of the total amplitude are not shown. [Bottom] The 
energies of the two localized modes compared with the local 
vibrational DOS of the outermost carbon atoms on infinite 
graphene edges. For H-passivated edges additional 1-D bands 
are found at ~380 meV (not shown). 
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FIG. 2: (Color online) [Top] The ZAZZZ structure. [Bottom] 
The three ADMs, shown in purple (c), green (d) and orange 
(e), respectively. Atomic displacements of less than 5% of the 
total amplitude are not shown. 



inserted. These two structures are chosen to mimic the 
experimental situation, and to test the generality of the 
trends found in our calculations. 



Using density functional theory [13l . Il4j | we shall 
demonstrate that strong local Joule heating occurs in 
systems of this kind for voltage biases and currents of 
the same order of magnitude as in the experiment by Jia 
et al. 12]. In order to have significant Joule heating two 
conditions are of importance. Firstly, localized vibra- 
tions are necessary in order to spatially concentrate the 
energy. Secondly, the electronic subsystem must couple 
strongly with the local vibrations and be locally out of 
equilibrium in order to provide energy to the vibrations. 
As we shall show, the structures shown in Figs. Q] and [2] 
indeed exhibit localized vibrational modes. These modes 
originate from "armchair-dimers", defined as C-C dimers 
coordinated in the same way as the outermost atoms of 
a non-passivated armchair edge. In what follows, we call 
these modes "armchair dimer modes" (ADMs). Next, we 
show that zigzag-armchair junctions exhibit strong local 
scattering of the current carrying electrons such that we 
can expect a local voltage-drop across the junction. Fi- 
nally, we estimate the heating of the ADMs in the two 
model systems and argue that the heating of the ADMs 
is the likely cause of the reconstruction of the edges ob- 
served in the experiment by Jia et al. [l2T |. 



Let us now address quantitatively the existence of lo- 
calized vibrational modes in systems with mixed zigzag 
and armchair edges. Using finite displacement calcula- 
tions [HI we find two such modes (see Fig. [J) for the ZAZ 
system with a vibrational energy of ~250 meV. These 
modes are strongly localized to the outermost atoms of 
the armchair edge: already for nearest neighbors the 
mode amplitude has dropped by ~85 %. The modes 
are truncated versions of the vibrational edge states that 
give rise to the quasi-lD band in the density of states 
(DOS) [16] of the infinite armchair edge (full black curve 
in Fig. [TJ see also Ref. Il7h . These two ADMs are en- 
ergetically localized since both the infinite zigzag edges 
(full grey curve in Fig. [T]) and bulk graphene (band edge 
shown as a dotted line) have a vanishing DOS above 
~200 meV and thus cannot cause energy broadening 
above this energy in the harmonic approximation. The 
modes outside the bulk band are only broadened by in- 
teractions with electrons and, at high temperature, by 
anharmonic interactions. 

For the ZAZZZ system we identify three localized 
modes (Fig. [2]) which again involve two-coordinated arm- 
chair dimers. For this case the vibration amplitude is 
significant not only on the armchair edges but also at the 
240° ZZ edge. 

Hydrogen passivation is predicted to play a role for 
the edge structure[l3|. However, we expect the ADMs 
to be localized even if they are adjacent to hydrogen- 
passivated edges. We find that the passivated zig-zag 
and armchair edges do not have vibrations with energies 
^250 meV in their DOS (dashed black and grey curves 
in Fig.[H respectively), in agreement with the empirical- 
potential calculation in Ref. and thus cannot pro- 
vide harmonic damping. We constrain our studies to 
non-passivated edge interfaces since there were no signs 
of hydrogen playing a role in the experiments by Jia et 
al. hj and we also expect hydrogen to desorb at much 
lower temperatures than where the breaking of carbon- 
carbon bonds takes place. The barrier for hydrogen des- 
orption from an armchair edge is lower compared to that 
of the outermost C-C unit, 4 eV (Ref. [13|) and 6.7 eV 
(Ref. [13]) > respectively. 

In addition to the existence of localized vibrational 
modes, the occurrence of local Joule heating also requires 
local scattering of electrons. By NEGF-DFT transport 
calculations 14. 21. 22| , we find that the scattering states 
localized at the zigzag edge are interrupted at the zigzag- 
armchair junction, illustrated by the square modulus of 
the transmitting electronic scattering states (eigenchan- 
nels [23j]) at the Fermi level, shown for the minority spin 
in Fig. [3J Ref. [24| shows, by extensive tight-binding 
calculations, that scattering can in general be expected 
at armchair-zigzag edge interfaces and that scattering in- 
creases with the length of the armchair edge. Since the 
same behavior is predicted for different systems with two 
different methods (ab initio vs. tight-binding), it seems 
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FIG. 3: (Color online) The absolute square of the minority 
spin left-to-right scattering states at the Fermi level for a large 
ZAZ system, corresponding to a transmission of 0.15. 



that back-scattering at armchair-zigzag interface indeed 
is a generic feature of these systems. 

We next estimate the heating of the ADMs in the 
two model systems. If the anharmonic couplings are 
neglected, the mean steady-state occupation, (n x )(V), 
can be calculated from the ratio of the current-induced 
phonon emission rate, 7 A m (V^), an d the effective phonon 
damping rate, 7damp(^0 (here V is the voltage drop 
across the scattering region). Since the investigated 
modes lie outside the vibrational bulk band, the damping 
due the bulk phonon reservoir vanishes. Assuming zero 
electronic temperature and energy-independent scatter- 
ing states within the bias window, the emission rate is 
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Here, tuo\ is the vibrational energy of mode A, M is the 
coupling of the mode to the electronic degrees of freedom 
calculated by finite difference techniques, and A^/^ is the 
electronic spectral density of left /right moving electrons, 
evaluated at the Fermi level by NEGF-DFT transport 
calculations 0, ED, El . The rate 7d amp (^) is 



Td amp 

(V) 



Tr[M A AM A A] , 



(2) 



where A = Al + Ar is the total electronic spectral den- 
sity at the Fermi level. The occupation (n A )(V) becomes 



where s A = 2T ^^^am>-a] LI ls a dimensionless heat- 
ing parameter that can vary from (no heating) to 1 
(maximal heating). By assuming that n x (V) is Bose dis- 
tributed, one can extract an effective temperature of the 
mode, T^ e (V), which we plot in Fig. |4] for the localized 
modes identified in Fig. [Hand Fig. [51 The calculated elec- 
tronic damping rates and heating parameters are shown 
in Table [B 
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FIG. 4: (Color online) The effective temperature as a function 
of bias for the localized modes for the ZAZ (full red and blue 
curves) and the ZAZZZ system (dashed purple, green and 
orange curves). The full black horizontal line indicates the 
uniform temperature where the decay rate of the outer C-C 
dimer of an armchair reaches 1 Hz. 



Mode 


hu>x (meV) 


FnLmp(tJ-eV) 


s x 


ZAZ, a 


248 


99 


0.77 


ZAZ, b 


239 


63 


0.29 


ZAZZZ, c 


244 


19 


0.63 


ZAZZZ, d 


238 


26 


0.38 


ZAZZZ, e 


234 


64 


0.28 



TABLE I: Mode characteristics. 



The mode temperatures should be compared to the 
uniform temperature, Td, needed to destabilize the sys- 
tem on a time-scale relevant to the experimental condi- 
tions, which we judge to be of the order of seconds. Thus 
we consider a corresponding rate of desorption of carbon 
dimers, q ~ 1 Hz. We estimate Td ~2500 K (see Fig. [4]) 
using the Arrhenius equation, q = v exp[—E a /(kBTd)], a 
characteristic attempt frequency hv = 100 meV (a typical 
phonon frequency), and an activation energy E a = 6.7 eV 

We can compare the heating of the localized ADM to 
the modes inside the bulk band by calculating their har- 
monic damping due to their coupling to the bulk phonons 
[28j |. and adding this to the damping by the electrons. We 
find that this damping is 1 — 100 times the damping due to 
the electronic couplings leading to temperatures typically 
below 1000 K even for a bias of 1 V. Since this tempera- 
ture yields desorption rates much lower than those seen 
in the experiments we conclude that the Joule heating 
of the harmonically damped modes cannot account for 
the reconstruction. Instead, as Fig. [4] shows, the ADMs 
reach a high temperature at much lower biases, and can 
thus provide a channel for local desorption. Furthermore, 
these modes are also more likely to be involved in the des- 
orption since they directly involve the desorbing dimers. 

A quantitative comparison with experiments would re- 
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quire a much more sophisticated theory, beyond the scope 
of this Letter. For example, one should consider the an- 
harmonic coupling of ADMs, and evaluate the nonequi- 
librium electronic distribution function, from which the 
actual potential profile in the sample can be extracted. 

Finally, let us investigate the role of different types 
of edge interfaces. Fig. Q] shows two examples: a 150° 
zigzag-armchair interface, and a 210° armchair-zigzag in- 
terface. Intuitively one would expect the 2-coordinated 
dimer directly at the 210°-corner to make this interface 
especially prone to reconstructions. This is confirmed by 
our calculations of their heating: we find that the 210° 
modes (modes (a),(c) and (d)) exhibit markedly stronger 
heating than modes associated with the 150° corner, the 
reason being their stronger coupling to the current car- 
rying electrons (s A closer to 1 for these modes). This sce- 
nario is consistent with two experimental observations by 
Jia et al.\v^. Firstly, they observe that a 150° interface 
survives even after massive reconstruction has occurred 
(Ref. [12 [. Fig. 2(D)). Secondly, certain armchair edges 



evaporate while others grow longer. The theory outlined 
here predicts that the evaporating armchair edges are the 
ones bordered by at least one 210° junction. Conversely, 
zigzag edges and armchair edges bordered by 150° junc- 
tions would be the stable edges, and would grow as the 
unstable armchair edges evaporate. 

In conclusion, we have demonstrated how specific C-C 
dimers can play a fundamental role in current-induced re- 
construction in graphene systems with mixed edges. We 
show that these dimers give rise to spatially and energet- 
ically localized modes, which give a natural explanation 
for the low onset bias for reconstruction observed in the 
experiments [l2j . Identifying the modes that heat up also 



allows us to make predictions of the overall behavior of a 
graphene sample under the influence of a current. Specif- 
ically, we predict that zigzag-armchair junctions with an 
angle of 150° would be more stable than the junctions 
with a 210° angle. We believe reasoning along these lines 
could contribute towards a quantitative understanding of 
other intriguing edge structures, e.g. the zigzag-reczag 
discussed recently in Ref. 
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